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We consider a gravitating spherically symmetric nonrelativistic configuration consisting of a mass- 
less chameleon scalar field nonminimally coupled to a perfect isothermal fluid. The object of this 
paper is to show the influence of the chameleon scalar field on the structure and evolution of an 
isothermal sphere. For this system we find static, singular and regular solutions depending on the 
form of the coupling function. A preliminary stability analysis indicates that both stable and un- 
stable solutions exist. For unstable configurations, by choosing the special form of the coupling 
function, we consider the problem of the gravitational collapse by applying the similarity method. 
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I. INTRODUCTION 

Over the course of the past decades, there has been substantial progress in cosmology and astrophysics supported 
by a growing quantity and quality of observational data requiring its theoretical description. Difficulties associated 
with such description are obvious: in considering various objects and processes in the Universe, one must deal with 
scales which differ by orders of magnitude. It forces the use of only approximate models and methods permitting 
considerable simplification in describing the structure and evolution of the Universe. In this respect, various scalar 
fields are ones of the most called-for objects which are being used in investigations of the Universe on all scales. They 
are relevant in describing the early inflationary Universe the present accelerated Universe 0, and the processes 
on scales of order of galaxies and separate stars. In particular, this has come to refer to attempts of describing the 
dark matter in galaxies [3] and to a consideration of compact configurations - boson stars supported by various scalar 
fields 0- 

In modeling boson stars, it usually assumes that they consist of complex or real scalar fields in their own gravitational 
field. However, it is also possible to imagine a situation when besides a scalar field there can exist other forms of 
matter in a system. This can be fermion fields Q, an electromagnetic field or ordinary matter. In the latter 
case, ordinary matter can interact with a scalar field either only via gravity or by use of nonminimal coupling as 
well. In the present paper we consider a case of the presence of such nonminimal coupling in a compact gravitating 
configuration. As applied to cosmology, the idea of nonminimal coupling between a scalar field and ordinary matter 
in the form of a perfect fluid was suggested in Refs. where the nonminimal coupling was used in describing the 

accelerated expansion of the present Universe. In this case the effective mass of the scalar field can change depending 
on the background environment. Because of it, such a scalar field was called a "chameleon scalar field." The idea 
that the properties of a scalar field could be influenced by the environment / matter surrounding the scalar field was 
studied earlier, in particular, in the papers 0, [l(| where the interaction between matter and the scalar field was used 
to model a dependence of fundamental coupling constants on the local environment. Further development of this idea 
can be found in the papers [TTI . fl2| , where the authors describe different forms of cosmological evolution within the 
framework of the chameleon cosmology. 

In progressing from cosmological to astrophysical scales, in a recent work [l3[ we considered the model of chameleon 
stars consisting of a scalar field nonminimally coupled to ordinary matter in the form of a perfect polytropic fluid 
having the equation of state p ~ p 1 , where p and p are the density and the pressure of the fluid, and 7 is some constant. 
For this system we found static, regular, asymptotically flat solutions for both relativistic and nonrelativistic cases 
for 7 = 2 and 5/3. For these values of 7, it was shown that the presence of the nonminimal interaction leads to 
substantial changes both in the radial matter distribution of the star and in the star's total mass. 

In the present paper we continue studying compact astrophysical objects consisting of a scalar field nonminimally 
coupled to the polytropic fluid having 7 = 1. This case corresponds to the isothermal gas sphere whose temperature 
is assumed to be constant along the radius. In the nonrelativistic case, there are known static singular and regular 
solutions describing such equilibrium configurations in weak gravitational fields [T3 |. Depending on the parameters 
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of the model, they can be gravitationally stable or unstable. In the latter case, instability leads to motion of matter 
that can be realized in the form of collapse. The gravitational collapse of an isothermal gas cloud is a well-studied 
problem. In different variants, this problem has been investigated frequently. In particular, there is a number of works 
devoted to the consideration of self-similar motions of an isothermal gas. We indicate some of them here. 

In considering the collapse, there has been a great deal of attention paid to the specification of initial conditions. 
Larson [l5[ and Penston |16l ] chose zero initial velocity and uniform density of the fluid. Such conditions provide a 
homologous inflow that evolves to a central region. This solution assumes that there exists the flow at large radii, 
or equivalently as t — > 0, which is directed inward at 3.3 times the sound speed and the density in the envelope is 
proportional to r~ 2 . It gives the density be 4.4 times the value appropriate to hydrostatic equilibrium. Shu [l7j 
criticizes this solution and suggests new self-similar solutions, which did not have these problems. He assumes the 
initial density profile in the form p ~ r~ 2 . Such a sphere is quasistatic, so that the infall velocities are negligible at 
the moment of protostar formation. The resulting initial configuration is the singular isothermal sphere, which is an 
unstable hydrostatic equilibrium. Shu showed that in this case the cloud collapse begins at the center, and the infall 
spreads outward. He called this solution the expansion-wave solution. Hunter [l8| has been continuing the work in 
that direction, and found a new class of similarity solutions which includes Shu's expansion- wave solution as one limit. 
He considered different initially unstable spheres and studied whether their collapses are well described by similarity 
solutions. 

In this paper we consider the model of a nonrelativistic configuration consisting of a perfect isothermal fluid in 
the presence of a massless chameleon scalar field. Our intention here is to clarify the influence of the chameleon 
scalar field on such a widely studied system as the gravitating isothermal fluid. The paper is organized as follows: In 
Sec. [H] we derive the equations for a static configuration with the arbitrary coupling function / (see the Lagrangian 
(fTJ) below). Using these equations, in next three subsections we consider various static configurations for different 
choices of /. For such configurations, we find analytical and numerical solutions, and also address the issue of their 
stability. Using as initial conditions unstable static singular spheres from Sec. Ill Al in Sec. IIIII we consider similarity 
solutions describing the collapse of such configurations. Our conclusions are summarized in the final section, where 
we also suggest some possible lines of further investigation. 



II. STATIC CONFIGURATIONS 



As discussed in the Introduction we consider a gravitating system of a real scalar field coupled to a perfect fluid. 
In the general case of strong gravitational fields, the Lagrangian for this system is 



L= — 



16ttG 



R 



V(tp) + f(<p)L T 



(1) 



Here (p is the real scalar field with the potential V(tp); L m is the Lagrangian of the perfect isotropic fluid i.e. a fluid 
with only one radial pressure; f(<p) is some function describing the nonminimal interaction between the fluid and the 
scalar field. An interaction similar to this was used to describe the evolution of dark energy within the framework of 
chameleon cosmologies [IJ EH • The case / = 1 corresponds to the absence of the nonminimal coupling. In this case 
the two sources are still coupled via gravity. 

The above Lagrangian was used by us in [l3j to construct static relativistic and nonrelativistic co nfig urations. In 



that case, the Lagrangian for the isentropic perfect fluid was chosen to have the form L m = p [li 



Lagrangian, the corresponding energy- momentum tensor is (details are given in Appendix in Rcf. [13j) 



20]. Using this 



= f [( P + p/c 2 )c 2 u lU k - 5'p] + d iV >d k <p - a* 



-d^ip - v(p) 



(2) 



where p and p are the mass density and the pressure of the fluid, u l is the four-velocity. Henceforth we will use the 
Lagrangian ([1} and the energy-momentum tensor @ both when considering static configurations and in the study 
of the collapse. It will be shown in Appendix lAl that the process of the collapse permits self-similar motions only in 
the case of a massless chameleon scalar field, i.e. when V(<p) — 0. Bearing this in mind, we restrict ourselves to the 
consideration of configurations supported by the massless chameleon scalar field only. 

Now we derive the equations describing static configurations. In generalized curvilinear coordinates the static 
spherically symmetric metric in the nonrelativistic limit can be written in the form (2~H 



ds 2 = e^c 2 dt 2 -dr 2 -r 2 <m 2 , 



(3) 
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where dVt 2 is the metric on the unit 2-sphere. Using this metric, the r-component of the covariant conservation law 
T" — (which is the only nonzero component in the static case) is 

dp 1 2 dv 

Tr = ~2 PC Tr- (4) 

In the nonrelativistic limit, in turn, the metric function v can be rewritten in terms of the Newtonian gravitational 
potential ip in the following form [21j : 

e" = l+?- (5) 

The corresponding equation for ip (the Poisson equation) can be found from the Einstein equations by using the 
transition to the nonrelativistic limit (2lj . Then, using the above expression for v and the energy-momentum tensor 
([2]) , and leaving terms of order 1 /c 2 only, we get 

1 d ( 2 dtp 
r 2 dr \ dr 



^Jr 2 ^)=^Gf P . (6) 



Next, introducing the mass 

dM A f 2 

— = iirfr p, 
dr 

one can find from the Poisson equation ©: 

d%jj _ GM(r) 
dr r 2 

Using this expression and taking into account ([5]), Eq. (0| takes the form 



i(- d l)=-**CfrV P) 



dr \ p dr 

This equation differs from the usual equation of hydrostatic equilibrium by the factor / = f(ip) reflecting the presence 
in the system of the nonminimal coupling. Correspondingly, a full description of the configuration under consideration 
requires an additional equation for the scalar field ip. In general form, it can be obtained by varying the Lagrangian 
((T|) with respect to <p as follows: 



_L d_ 

i—g dx % 



dx> 



-99 a u 



L ^ 

d<p 



where the metric gik is taken from ([3]) . This field equation with the above Lagrangian for the perfect fluid L rn = p, 
and the metric ([3]) with ([5]), gives the following scalar field equation: 



ji , 2 .j yf 



ip" + -iff = -p-f- , (8) 
r dip 

where the prime denotes differentiation with respect to r. Eqs. ([7]) and ([8]) should be supplemented by an equation 
of state. In the case of an isothermal fluid being considered here, we have 

p = Kp (9) 

with K — a 2 , where a is the speed of sound in the fluid. Using this equation of state and introducing the dimensionless 
variables 

e " = r e = z , <t> = ^^ (10) 

where L = ^K/{AivG\) has dimensions of length, A is an arbitrary constant, Eqs. (J7J and © can be rewritten as 
follows: 

I d ( t2 dr\ 



¥di{ e Tz)=^ ^ 
edt v dc. 
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In the absence of the scalar field, i.e. when / = 1, only the first equation remains which has the following solutions fl4j : 

(i) an analytical singular solution in the form e~ r ' = 2/£ 2 which diverges at the center of the configuration as £ — > 0; 

(ii) numerical regular solutions which satisfy the boundary conditions 77 = 0, drj/d^ — at £ = 0. These solutions 
approach asymptotically, as £ — > 00, to the singular solution from (i). Both sets of solutions imply that e~ v — > only 
asymptotically as £ — > 00. That is, such solutions describe only infinite size configurations. 

It will be shown below that the presence of the chameleon scalar field permits the existence not only singular and 
regular solutions similar to those described above, but also new finite size, regular solutions. If we are going to seek 
such solutions, it is obvious that a crucial role will be played by the form of the coupling function /. In the next three 
subsections we consider several types of solutions. In doing so, we will track the behavior of the energy density Tq 
from In the nonrelativistic approximation using the above dimensionless variables (|10[) . the total mass density 
has the form 



^T°/c 2 = A 



1 2c 2 \d£ 



(13) 



The second term on the right-hand side of this expression is proportional to the square of the ratio of the speeds of 
sound and light, (a/c) 2 . It is obvious that this term is small compared with the first term everywhere just except 
perhaps the vicinity of the points where both terms tend to zero simultaneously. 

Below, we demonstrate a few examples of solutions with different choices of the coupling function /. In making the 
choice of /, two approaches are possible: (i) One assumes / to be an arbitrary function of the scalar field (j> whose 
form is chosen to satisfy the requirement of obtaining the solutions needed. As an example, in Sec. Ill Al one of the 
simplest choice of / in the form of the exponential function, / = e~*, is presented. In this case, it is possible to get 
analytical singular and numerical regular solutions, (ii) One can suggest / to be a function of the radial coordinate, 
/ = /(£). As an example, we choose below power functions for / giving either infinite size singular solutions (see 
Sec. Ill Bp or finite size, regular solutions (see Sec. Ill Cp . The solutions for the scalar field cf> = </>(£) obtained in this 
case define the parametric dependence / = /(</>(£))■ 



A. Singular and regular solutions with / = e 

In this section, we seek a solution of the system (|TT]) - (fT^j) in the particular case when the coupling function / is 
chosen in the form 

/ = foe-*. (14) 

The parameter /o > can be absorbed by introducing the rescaling \/7o£ — > £ that allows to put /o = 1 in further 
calculations. But we will bear in mind that /o can be always restored in final expressions by using the above rescaling. 
The above choice of / allows to find the following analytical solutions 

«- = *_ h (At)-f W 

where A, B are integration constants. These solutions describe a configuration which is singular at the origin of 
coordinates, where, as £ — > 0, the mass density from (fT3"]) diverges as l/£ 2 , and regular asymptotically, as £ — » 00, 
where p t — > 0. 

The coupling function / from (|14[) permits also the existence of regular solutions which can be found numerically. 
In this case we can choose the parameter A from (|10p to be the central density of the fluid, A = p c . With this 
normalization the system (|li p -(|12 jl takes the form 

]_d_ ( pdn 

e V= dt 



72 17 U 37 = ( 16 ) 



ilf^|=e-H (17) 



We will look for a solution of this system which satisfy the following boundary conditions at £ = 0: 

77 = 0, ^ = 0, = (^ o = const, f§ = 0. (18) 
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Taking these conditions into account, one can show that solutions of Eqs. (|16|) and (fTTj) are related through . 
This allows the possibility of reducing the system p^|) - ([T7| to one equation: 



1 d 



e \^ 



2 d V 



-2r/-4>o 



which, after the rescalings £ = £ e^°/ 2 / 's/2, r\ = 9/2, takes the form 



1 d 



C 2 d( [ C dC 



,d6 



(19) 



This equation is an analogue of the Lane-Emden equation for the isothermal fluid without a scalar field [14j |. Thus, 
Eq. ([T9| [or equivalently Eqs. (JT6J) and (fT7|) ] together with the boundary conditions fp~8|) describes the regular isother- 
mal gas sphere in the presence of the chameleon scalar field. Without the scalar field, such configuration has the 
central density p c . In the presence of the scalar field, it follows from Eq. (fl3|) that the central density of the configura- 
tion under consideration is now defined by the value of the function / at £ = as well. In the case when the coupling 
function has the form (|14p . we have from (|13[) : pt c — p c & ■ It implies that the configuration under consideration 
will have a greater or smaller concentration of matter at the center depending on the sign of 4>o-i and correspondingly 
different distributions of the mass density along the radius. The size of such configuration, as in the case without the 
scalar field, will also be infinite. The corresponding asymptotic solutions of the system (|T5|) - (|T7)) are: 



V = 



ln£ 



Pi 



<T as £^ 



As in the case of an isothermal fluid without a scalar field, here we will consider a configuration embedded in an 
external medium of pressure p e . For such static configurations without a scalar field, in the works of Ebert [22j and 
Bonnor [23j], the issue of their stability was investigated. Following these works and taking into account the equation 
of state from (191. we represent the central density of the fluid p c of the configuration under consideration as follows 
(see also Ref. [18() : 



Pc 



p e 



(20) 



Then, taking into account the total pressure pt, which includes both the pressure of the fluid and the contribution 
from the scalar field, and follows from ©, p t = fp+ (d(p/dr) 2 /2, in terms of the dimensionless variables (TIT)]) , we have 



Pi 



a 2 p c 



di 



Substituting this expression in (|20[) and taking into account that at the outer boundary of the configuration, where 
£ = £ e , the total pressure is equal to the external pressure, i.e. pt{^e) — Pe, we finally have 



Pe 

P c = ~~2 



Then one can find that the total mass of the configuration is 



(21) 



M 



G^Gpe) 1 / 2 



1/2 



The maximum possible M , given p e , or the maximum possible p e , given M, follows from the fact that the function 
£ 2 \J j [(/>(£)] e _? '(5) + l/2<p' 2 dril 'd£ attains a maximum value at some £ = £ mQX . In the case when / is chosen in the 
form (|14p . the maximum value of £ is £ max = 5.235 e^ 2 . One can see that this value depends on the initial value of 
the field 4>o, an d differs from the case without a scalar field when ^ ma x = 6.451 [3, [23|]. The maximum value of the 
radial coordinate, £, ma x — 5.235 e^°/ 2 , corresponds to the mass of the cloud having the critical value 



M crit = 0.52 



G 3 / 



:P< 



-1/2 



(22) 



This value is factor of 0.44 smaller than the critical Bonnor-Ebert sphere. Notice that the critical mass does not 
depend on the initial value of the scalar field 4>q and is determined just by the external pressure p e , as in the case of 
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FIG. 1: Total pressure distributions of bounded isothermal spheres in the presence of the chameleon scalar field with the coupling 
function / taken from (|14|l . All nonsingular curves are plotted at the initial value of 4>o = 0. The curve marked "critical" 
corresponds to the configuration with the critical mass given by Eq. (|22[l . The spheres which are less centrally concentrated that 
the critical configuration are gravitationally stable; those which are more centrally concentrated are gravitationally unstable. 
In the limit of infinite central concentration, the latter spheres approach the singular solution. 



the critical Bonnor-Ebert sphere. The relation between the central density and the total density at the outer boundary 
of the configuration, p c j p t — f^ 1 e v , corresponding to this critical mass, can be written as p c J pt(£max) — 19.724 e^°. 
Then for clouds with masses less than M cl a two types of configurations are possible: a larger stable one with a range 
of density less than 19.724 e^°, and a smaller unstable one with a range of density greater than 19.724 e^°. For both 
cases £ e < ^max- Configurations with £ e > £, max are always unstable, and with £ e = ^ mQX are marginally stable. 

In the limit of configurations with infinite central density described by the singular solutions (|15p . the unstable 
equilibria approach the singular sphere which has the following density and mass distributions: 

*(r) = ^|r->, M(r) = jf^r. (23) 

This singular solution, truncated at a boundary pressure p e , has a total radius R and a total mass M(R) given by 

1/2 

which are, respectively, factors of 0.58/ and 0.24 differ from the critical Bonnor-Ebert sphere. Here we restored 
the parameter /o from (I14[) whose value allows varying the radius of the configuration for the given mass. The plots 
of total pressure distributions of several configurations described above are shown in Fig. [T] 



B. Singular solutions with / = /(£) 

In this and the next subsections, we seek solutions of the system (fTT|) - (fT5|) starting not from a straightforward choice 
of the specific form of the coupling function /, but by choosing an appropriate form of the product /e - '', defining the 
total density (fT3|) . In paper il3j, we already found the analytical solutions for a massless scalar field in a particular 
nonrelativistic case of an incompressible fluid. Below we will obtain analytical solutions for the case of the isothermal 
fluid being considered here. 



7 



Choose first the above product in the form 



where B,p are arbitrary positive constants. This choice gives the positive-definite total mass density p t from JTe 
which is singular at the origin of coordinates. The choice (1231) allows to find the general solution of Eq. (fTTj) in the 
form 



ft-" = §■ ( 25 > 



* = c ' + T + 0-X-») ' 



In the case when p = 2, 3, we have 



m = Ci + y + Sln^ % = C?i + y - I (1 + InO . 

Here C±,C2 are integration constants. The scalar field equation (fT2j) . in turn, can also be integrated analytically 
giving the expression for 0' 2 . 

Of special interest is the case of p = 2. In this case one can get a solution similar to the singular solution without a 
scalar field. Indeed, choosing the integration constants as C2 = and C\ = In (1/2), we obtain the following solutions: 



i.'V ^ = ^r^ h ^" 1*.. 



In the case when B = 2we return to the singular solution without a scalar field from (i) [see after Eq. ([T2])]. The values 
of the constant B ^ 2 imply that the scalar field changes along the radius that provides changes in the distribution 
of the mass density from (IT3| . 

It can be also shown that, putting B = l,p = 2 in (|2"5)l , solutions (f2l)|) correspond to the singular solutions from 
Sec. Ill Al if one chooses the parameters from expressions (fT5)) to be A = 2, B = 0. Thus we have 

It is obvious from two first expressions that the parametric dependence / = /(</>(£)) gives the expression / = 
coinciding with the coupling function from (|14p. In Sec. IIIII below, a self-similar motion of matter with such choice of 
/ will be considered. 

C. Regular solutions with / = /(£) 

Bearing in mind that we seek finite size, regular solutions, let us try to choose such an expression for fe~ v that 
provides the required solutions. In doing so, we will look for such solutions of Eqs. (fTTj) and (fT2|) which have a value 
of the total mass density pt from (fT3|) at the center of the configuration, at £ = 0, equal to the central density of the 
fluid, A = p c , i.e. p t o = p c , and at the outer boundary of a cloud, at £ = £1, we require that pn — 0. One of the 
simplest variants is to take the following power-law dependence 

/e-" = l + a^, (27) 

where a,/3 are arbitrary constants. Taking all the above into account, we choose a to be negative, and fi > 0. 
Substituting expression (|2T|) in Eqs. (fTTj) and (fT2"j) . we find their solutions in the form 

V = li 2 + ta 1 ° , ^ 2+ff , (28) 



6" (0 + 2)08 + 3) 



( 2 = _i!Ji + _^ 
dey 9 \ ^+3 S 



1 , W + 3 ) 1 3q t 
/3 + 4 2(/3 + 3) ? 



(29) 



Next, from expression (|27|) . let us determine the point £ = £1 in which /e^ = 0: £1 = (— a) -1 /' 3 . Requiring that in 
this point the derivative 4>' also be equal to zero, we substitute this value of £1 in Eq. (121)1) and find the corresponding 
value of a: 

f /3Q8 + 4) 
1 18(0 + 3)2 
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Using this expression, we get 

( /?(/? + 4) ^ ' 2 
U8(/3 + 3) 2 
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This point can be interpreted as the outer boundary of the configuration, where the total density from Eq. (|13[) is 
equal to zero. In this case the spheres described by Eqs. (fTTj) and (fT2"j) are a one-parameter family, with the size of 
£i depending only on one parameter j3. As (3 — > oo, the value £i — > 3-\/2, and as /3 — >■ 0, we have £i ~ l/V? — * 00 ■ 
Thus, by changing the value of the parameter /?, we can change the size of the configuration in a wide range at the 
given p c and K. 

Note that the scalar field energy density which is proportional to 4>' 2 from (|29l) , can take both positive and negative 
values along the radius depending on the value of the parameter (3. This is related to the structure of the scalar 
field equation ((SJ) which, at positive pressure p and a specific choice of /, can effectively correspond to a scalar field 
equation for a ghost scalar field. Since the form of / = f is defined parametrically through the solution for the 

scalar field <j> = </>(£), and depends on the value of the parameter /?, then a situation may occur where the scalar field 
equation (J5J demonstrates the ghostlike behavior with <f>' 2 < 0. However, since in the nonrelativistic limit the main 
contributions to the energy density (|13p come from the first term on the right-hand side, which is always positive, 
then the total energy density remains positive as well. 

Returning to the dimensional variables from (|10[) . the radius of the configuration under consideration is given by 

1/2 



In Sec. IIII1 where we will consider the process of the collapse of the singular sphere from Sec. Ill A[ the speed of sound 
is chosen to be a = 0.2 km s _1 . If we use the same value in the estimation of the size of the regular configuration 
which is considered here, we have 

R = 2.18 x 10 7 — cm, 
yfp~c 

where p c is the central density of the configuration in g cm -3 . 

It is seen from the obtained results that the isothermal fluid in the presence of the nonminimally coupled scalar 
field with the coupling function / defined through Eq. ([27]) permits the existence of finite size, regular solutions. This 
situation differs from the classical problem without a scalar field when infinite size solutions are the only possible ones 

B n 

Let us now consider the issue of the stability of the solutions obtained. Following [24], using Eq. (fTTj) . we introduce 

a quantity 



Because of the nonrelativistic character of the problem, we neglect the influence of the scalar field. Next, using 
expression (|13[) , we define the average density as 

M M 

Pt 



(4/3)ttR 3 (4/3)ttL 3 cJ 3 
For the density concentration, we have 



By combining two last expressions, we find 



Pt 3/Lii 

M 



Taking into account the expression for L from (fTO)) . we obtain 

Pc = 



One can see from this equation that p c drops with increasing M, which is the abnormal behavior for a stable star in 
equilibrium. Thus, regardless of the value of the parameter /?, the choice of the coupling function / in the form (|27|) 
gives only unstable configurations. 
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III. SELF-SIMILAR COLLAPSE OF THE SINGULAR SPHERE 



In this section we consider the spherical gravitational collapse of the singular configuration investigated in Sec. Ill Al 
Such problems have been considered frequently for polytropic spheres both in relativistic and nonrelativistic cases. 
In doing so, two main approaches are being used: (i) One solves the system of hydrodynamic partial differential 
equations; (ii) A self-similar motion of polytropic matter is considered. 

Here we study self-similar motions of the matter of the singular configuration considered in Sec. Ill Al In the absence 
of the chameleon scalar field, such problem for an isothermal sphere have been considered by many authors (see the 
Introduction). We obtain below the similarity solutions describing the collapse of the singular sphere in the presence 
of the chameleon scalar field. These solutions will be compared to the known solutions found in Shu's paper (f7| . 

In Appendix [X] we derived the following set of nonrelativistic hydrodynamic equations which take into account the 
influence of a chameleon scalar field: 



d 2 tp 
dr 2 



dM 2 
+ Anfr pv 

dM 
dr 

dv dv 
dt dr 
2ckp dV_ 
r dr dtp 



0. 

4vr/r 2 p, 

I dp GM 
p dr r 2 
Of 



= -P 



dip' 



(31) 
(32) 

(33) 

(34) 



Here M = M(r, t) is the total mass inside radius r at time t, v the radial velocity, p the pressure, and p the density. If 
the coupling function is chosen in the form / = , then the corresponding system of self-similar Eqs. (|A.21[) - (|A.23|) 
from Appendix [Al can be written as follows: 



da a 

dx (x — u) 2 — 1 

du 1 



e" a- 2 V (x-u) 



dx (x — u) 2 — 1 
d 2 4> 2 dej) 
dx 2 x dx 



b a(x -u) h-r 

x dx 



dx 
(x - u), 



(x - u), 



ae 



(35) 
(36) 
(37) 



The similarity variables involved were obtained from expressions ()A.I4|) in the case when 7 = n = 1, k = K = a 2 . 
This gives 



— , v{r,t) = au(x), p(r,t) 
at 



a(x) 
4nGt 2 ' 



p(r,t) 



AirGt 2 



a(x), M(r,t) 



aH 



m{x), ip(r,t) = —j== 
G \/AttG 



(j>{x), 



(38) 



where t = defines the instant when the mass of the core, M(0,t), is zero. The instant t = corresponds to the 
instant of core formation for the collapse problem. In the case of collapse the variables t, x, m are positive while u is 
negative. 

In terms of the above similarity variables we have the following expression for the mass function m(x) from ()A. 151) : 



m(x) — e ^x 2 a(x — u). 



(39) 



Thus we have the system of three Eqs. (|35]) - (l37l) and relation ([39]) describing the self-similar collapse of the configu- 
ration under consideration. 

An exact analytic solution of Eqs. ([3"5"]) - (|3"Tl) is given by the static solutions (fT5|) in the particular case of B = 0: 



u = 0, a = — , (j> = In (Ax) , m 



(40) 



In dimensional units, these solutions correspond to the time- independent singular solution (|23[) with /q = 1. As in 
the case of the problem without a scalar field considered by Shu in [17| . in our case this singular solution is the only 
hydrostatic solution allowing self-similarity. Following [l7|, we will use this solution as the "initial state" when the 
process of collapse starts. 



10 



Below, we will seek numerical solutions of the system (1351) - (1371) . However, this numerical study needs some caution 
because the solution can pass through critical points, as well as in the case of the absence of a scalar field [17j. The 
existence of the critical points assumes that the denominator of Eqs. (|35|) and f[36]) vanishes: 



(x - u) 2 - 1 = 0. 

Substituting this expression in Eqs. (|35|) - ([37l) . one can obtain the following set of equations: 

-j-+e =0, 

d 2 (j) 2d<j> A V ' 

— — H — = ae~v. 

dx 2 x dx 

Solutions of these equations make the numerators of Eqs. (|35|) and (|36|) equal to zero, thereby providing the regularity 
of the solutions at the critical points. The particular solutions of the system (|41]> are 

a = 2e 2/x , (/)=- + 2lnx. 

x 

One can show that these solutions do not satisfy the original Eqs. (|3"T j) -(f34 |) . Henceforth we will only consider the 
solutions which do not pass through the critical points. 

To find a solution of Eqs. (J35J- (J3TJ) , it is necessary to specify the corresponding boundary conditions. When 
considering a problem of collapse, it is reasonable to assume that the fluid velocity is negligible at the "initial instant," 
i.e. that u — > as x — > oo. Then one can show that solutions of the system (f3"5]) - (f3"7) which have this property have 
the following asymptotic behavior: 

a^A/x, u ~ Uooe l ^ 2x2 \ cj> ~ In (Ax) , m~x as x -> oo, (42) 

where < is a value of the velocity u as x — > oo. The collapse begins from the initial inhomogeneous distribution 
of the matter in the form 

ftMJ-^r-. (43) 



which corresponds to the singular static sphere from (|23|) . This sphere being initially in unstable hydrostatic equilib- 
rium, can spontaneously generate inflow at every radii when ^ 0. Several solutions of the system (|35p - (|37[) with 
different initial values of the velocity are presented in Fig. [2] 

Shu starts his similarity solutions with the initial singular state of the type ()43|) using as initial conditions the 
asymptotic expressions for the velocity u ~ —(D — 2)/x and for the initial density distribution p(r, 0) ~ Dr~ 2 
depending on the value of the arbitrary parameter D. In the limiting case D — > 2, Shu returns to the singular 
sphere which is the solution of static equations. As a result, he obtains the limiting self-similar solutions, called the 
"expansion- wave collapse solutions." In our case we also have similar solution shown in Fig. [5] by the thick solid line. 
We will describe this solution in detail in Sec. IIII Al 

Let us now consider the behavior of collapsing solutions of the system (I551 - ([3T)) in the limit of x <C 1. Assuming 
here supersonic flow, u 2 3> 1, we obtain the following solutions near the origin: 



1 -b/x q+ \/q 2 + 2b b 
to ^ to , a ^—e b ' x , u -> -2 j <P^-~. (44) 



x 



Here 6, q are constants whose value is determined from the solution, and depends on Uoo , and too = q + \J q 2 + 2b 

is the value of the mass at the center. When 6 = 0, i.e. in the absence of the scalar field, we return to the known 
solutions from Ref. 1_7|: to — > too, a —> (mo/2x 3 ) 1 ^ 2 , u — >• — (2mo/x) 1 ^ 2 . 

Following [T3|, let us calculate the ratio of the gravitational acceleration to the pressure gradient acceleration for 
the fluid: 

GM a 2 dp to d In a too 
r z p or x z dx b 

which tends to the constant at finite t as r — > 0. This situation differs from Shu's result, where the above relation 
diverges that is quite expected behavior for the collapse problem. The obtained result is obviously related to the 
presence in the system of the scalar field whose behavior becomes dominant in the limit x — > 0. In this case one can 
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FIG. 2: Similarity solutions for the velocity u with different initial values of ttoo. The thick solid curve gives the collapse 
solution for the singular sphere (|43|l which is initially hydrostatic ( "expansion- wave collapse solution"). 

suppose that it is needed to perform a more correct study of the problem assuming the presence of the scalar field in 
the expression for the energy density TJ from (|A.2[) which is not taken into account in the present model because of 
its nonrelativistic character. 

The constant mo, which represents the reduced mass of the core, is expressed through the above constants b, q 
which, in turn, are defined implicitly by the constant Uoo from (|42j) . but not by A which, as it is obvious from the 
structure of Eqs. (I3"5j) - (l3"71) . does not influence the solutions. The values of b, q, and correspondingly mo can be found 
by numerical integration with the starting conditions (|42|) determined at large x. Numerical calculations show the 
following linear dependence which is approximately valid at {u^l < 10: mo ~ 0.491 + \uqo\. 



Following Shu |17j , in this section we consider similarity solutions starting from the initial static singular distribution 
of matter in the form (|43p . Then, letting the initial velocity Uoo from P2"]l equal to zero, the solutions obtained below 
can be interpreted as follows: we start from the initially unstable singular sphere (I43[) with zero initial velocity u 
at the instant t = 0. The collapse begins at the central region as a small perturbation of density spreading in the 
external envelope at the speed of sound a up to the outer boundary of the configuration, r — R, defined by Eq. (|24l) . 
Behind the front of this expansion wave, the motion of the matter is described by Eqs. (|35[) - (|5T|) whose solutions are 
sought in the range + < x < 1. The value x ~ 1 corresponds to the outer boundary of the configuration where 
r = R. The time required for initiation of the collapse over the whole configuration is T = R/a. During the collapse, 
in the neighborhood of the center of the configuration, the accumulation of mass from the initial value mo — 0.491 to 
m = 1 takes place. The latter value represents a total reduced mass contained inside x — 1 and equal to the original 
equilibrium value from pi))) . The physical mass M contained in the expansion wave and in the core increases linearly 
with time according to Eq. (I38[) . As well as in Shu's work, the mass contained in the core at any instant t always 
comprises about 49% of the total mass contained within r = at. 

Let us now illustrate the process of collapse described above by using numerical calculations with a specific choice 
of the value of the speed of sound, a — 0.2 km s _1 [llj ■ In Ref. [llj the external pressure was chosen to be p e /k = 
1.1 x 10 5 cm~ 3 K (here k is the Boltzmann constant) which defines the size and the total mass of the configuration 
as R — 1.6 x 10 17 cm ~ 0.052 pc and M — 0.96A/©, respectively. It allowed Shu to identify such an object as a Bok 
globule embedded in an H II region, and to use it as an initial state from which the process of collapse started. If we 
want to get the same value of the mass as in [l7j , then by comparing our expressions for the size of the configuration 



A. Expansion-wave collapse solution 
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logic ( r ) log, ( r ) 

FIG. 3: Expansion-wave collapse solution for a 0.96Mq singular sphere with a — 0.2 km s _1 and p e /k = 1.4 x 10 4 cm -3 K. The 
initial radius of the outer boundary is shown on the left panel by the vertical dashed line. The total density pt and the velocity 
— v profiles at t = 0.25, 0.50, 1.00, and 2.30 x 10 13 s are shown (from top to bottom on the left panel, and from left to right on 
the right panel, respectively). The dimensions of r,p t , and v are cm, gem -3 , and cm s -1 . 

and its mass, given by Eq. with the corresponding expressions from [l7|, given by Eq. (3), it can be seen that 

in our case one needs to take the pressure p e 8 times smaller than Shu's value. It will lead, in turn, to the growth of 
radius R which is now y/8 times larger than the value from . 

The collapse begins in the neighborhood of the center of the configuration at the instant t — 0. The expansion wave 
moving through the configuration with the speed of sound a, reaches the outer boundary at R = 4.5 x 10 17 cm ~ 
0.147 pc (that is about 1000 times the distance of Neptune from the Sun) at the instant T = 2.3 x 10 13 s ~ 7.29 x 
10 5 yr. The profiles of distributions of the density and the velocity in the neighborhood of the outer boundary of 
the configuration, r — R, are shown in Fig. [3] at different instants of time. Notice here that in performing numerical 
calculations, it is convenient to introduce new exponential function a = e~ a . 

By comparing the results presented in Figs. [5] and [3] with Shu's calculations shown in his Figs. 2 and 3, it can be 
seen that the presence of the chameleon scalar field with the coupling in the form (TT4")) does not result in qualitative 
changes of the picture of the collapse. Changes touch upon just some quantitative characteristics: starting from the 
task of obtaining the mass of the collapsed object of order a stellar mass, we have to take the external pressure 8 
times smaller, and the size of the initial configuration y/8 larger, compared to Shu's problem. Near the center, the 
process of collapse is analogous to the problem from Ref. [l?} when the total mass density pt ~ r~ 3 / 2 and the velocity 
u ~ r -1 / 2 [see expressions (|44j)] that provides the constant in time central accretion rate. 

IV. SUMMARY AND CONCLUSIONS 

In this paper, we have studied nonrelativistic gravitating configurations consisting of an isothermal fluid nonmin- 
imally coupled to a chameleon scalar field 4>, and described by the Lagrangian ([T]). The aim of this work was to 
clarify how the scalar field affects the distribution and motion of matter in such configurations. The crucial role here 
is played by the form of the coupling function / = /(</>). In considering these problems we restrict ourselves to the 
particular case of a massless chameleon scalar field. Besides allowing us to find some analytical solutions for static 
configurations, it allows us to consider the problem of the gravitational collapse by applying the similarity method. 

Namely, in Sec. Ill Al the static singular and regular solutions were found for the case of / = e~^. It was shown that 
there exist analytical singular solutions (|T3|) similar to known solutions for an isothermal sphere without a scalar field 
(13 ) . The study of the regular solutions indicates that the properties of the configuration under consideration depend 
substantially on the central value of the scalar field <t>^: (i) Depending on the sign of (f>o, the configurations will have 
a larger or smaller concentration of matter at the center; (ii) For a configuration embedded in an external medium 
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of pressure p e , the hydrostatic equilibrium has the critical value of mass (|22|) which is factor of 0.44 smaller than the 
critical Bonnor-Ebert sphere [22l . 23]; (iii) There is a critical size of the configuration, £ = £ mox , corresponding to the 
above critical mass, whose value is determined by the value of 4>o- In this case configurations with sizes less than 
£max may be either stable or unstable depending on the central density, and configurations with £ > £ maw are always 
unstable. 

Next, in Sees. Ill Bl and lll Cl the energy density of the system is chosen to be a power function of the radial coordinate 
£ in the form (|25p and (|27[) . For these cases, singular and regular solutions were found. The regular solutions permit 
the existence of finite configurations whose size is determined just by one free parameter - the exponent /? in expression 
(|27|) . This situation differs from that in problems without a chameleon scalar field when only solutions describing 
infinite size configurations do exist. A simple stability analysis indicates that, for the choice of / in the form (|2"T|) , the 
regular, finite size solutions are unstable. 

The above studies indicate that systems under consideration can be stable or unstable. In the latter case, a question 
arises as to the possible motion of matter in its own gravitational field. As an example of such motion, in Sec. IIIII 
the self-similar collapse of the singular sphere with function / chosen in the exponential form (fT4|) was considered. 
Without a scalar field, the similar problem was investigated by Shu in [13]. That is why we basically compared our 
calculations with Shu's results. It was shown that the qualitative behavior of the solutions remains the same as in 
the case without the scalar field: the process of collapse starts from the initially singular sphere (|23|) truncated at a 
boundary pressure p e and having the radius and the total mass given by expressions (|24p . The collapse is initiated by 
an initial density perturbation near the center that propagates outward to the edge of the configuration in the form 
of the expansion wave with the speed of sound a. The presence of the scalar field gives some quantitative differences 
which are obviously connected with the fact that the initial density distribution (I43|) differs from Shu's distribution 
by the factor 1/2. As a result, starting from the task of obtaining the mass of the collapsed object of order a stellar 
mass, it is necessary to set the external pressure 8 times smaller, and the size of the initial configuration y8 larger, 
compared to Shu's problem. 

As in Shu's problem, in this paper we consider the collapse of the singular isothermal sphere bounded by some 
artificial external pressure. Without the scalar field, inclusion of such external pressure is the only possibility of 
obtaining finite size configurations with finite masses (both for singular and regular cases). As shown in Sec. Ill CI 
in the presence of the chameleon scalar field interacting with the fluid in the form (|27[) . it becomes possible to find 
regular solutions describing more realistic finite size configurations without the external pressure. Such configurations 
being initially unstable, in principle, might be used as initial states when considering the collapse of isothermal clouds 
in the presence of the scalar field. 

Another possibility of obtaining finite size solutions can come from a consideration of Eqs. (fTTj) and (IT21 as a 
system which is equivalent to two interacting scalar fields rj and <p. In this case it is possible to introduce the effective 
Lagrangian of the system in the following form: 



where the effective potential V(rj, 4>) — — e ~ r '/(^)- By varying this Lagrangian with respect to 77 and </>, one can obtain 
Eqs. (fTT|) and (fT2]l . Then the problem of searching the finite size solutions to the field equations amounts to finding 
such a form of the coupling function / that provides the required solutions, if any. 



The author is grateful to the Research Group Linkage Programme of the Alexander von Humboldt Foundation for 
the support of this research. 



In this appendix we derive equations describing self-similar motion of a perfect isothermal fluid in the presence of 
a chameleon scalar field. To do this, we start with the energy-momentum tensor @. Let us choose the metric in the 
Newtonian approximation in the form 
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Appendix A: Derivation of the self-similar equations 




(A.l) 
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where dfl 2 is the metric on the unit 2-sphere, and ip is the Newtonian gravitational potential. Using this metric and 
taking into account the expressions for the components of the four-velocity 



u 



/goo 



Ul = 



9n 



v/c 



/500 a/1 - (v/c) 2 ' 



where v is the ordinary three-dimensional velocity, one can obtain the following components of the energy-momentum 
tensor @ in the nonrelativistic limit: 



1 



V(<p), 



7f =Ti = -fp+ -ip' 2 +V(p), ^ = -T° = fpvc. 
Next, using the law of conservation 

we get two equations: The i = component corresponding to the hydrodynamic equation of continuity 



(A.2) 
(A.3) 



The i = l component (Euler's equation): 

l flTf dTl 
c dt dr 



dt 



m-T°) 



1 9V 

c 2 dr 



•if) 



Substituting expressions (|A.2|) and (|A.3|) in these equations, we have, respectively, 



9p dp 
dt dr 



dv 
dr 



and 



\ 9t 9r p dr dr J 



-pv 



pdf 
p dr 



(A.4) 



(A.5) 



where the prime and the dot denote differentiation with respect to r and t, respectively. In the case when / = 1 and 
the scalar field is absent, Eq. (IA.4I) reduces to the usual equation of continuity, and Eq. (| A.5|) to Euler's equation for 
the spherically symmetric case. 

The equation for the scalar field ip coming from the Lagrangian ([T]) is 



1 d_ 

i—g dx % 



-99 



ik 



d<p 
dx k 



dV d£_ 

dip m dip 



Using this field equation with the Lagrangian for the perfect fluid L rn = p, and the metric (|A.1|) , we get the following 
scalar field equation 



<P 



-tp 
r 



dV 

dip 



-P— 
dp' 



(A, 



In obtaining this equation, we neglect the term ip/c 2 , in the spirit of the nonrelativistic character of the problem. In 
this approximation, Eq. (|A.6[) is static, i.e. one assumes that perturbations of the scalar field propagate with infinite 
velocity. Substituting Eq. (|A.6j) in (|A.5j) . we get the usual Euler's equation 



dv dv 1 dp dip 
dt dr p dr dr 



Next, introducing the mass M contained within a sphere of radius r 

4wfr 2 p, 



dM 
dr 



(A.7) 



(A.8) 
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one can rewrite the equation of continuity (IA.4I) in the following form: 

^ + 4nfr 2 pv = 0. 

Substituting expression (IA.8I) in the Poisson equation 

LL (r 2 ^\ =4nGf P 
r 2 dr \ dr J ' 

we obtain the usual equation for the Newtonian potential 

dtp _ GM 
dr r 2 

Substituting this in (| A.7|) . we get the final form of the set of equations describing the dynamics of the system under 
consideration: 

^+4TTfr 2 pv = 0, (A.9) 

^ = 47r/r 2 p, (A.10) 

dv dv 1 dp GM 

at or p or r z 

d ^t + l d JL- d -L = (A12) 

dr 2 r dr dip dip 

Here M = M(r,t) is the total mass inside radius r at time t. These equations are invariant under the time reversal 
operation: t — > —t, v —> —v that allows us considering only the range of < t < oo. Notice here that in the absence 
of the nonminimal interaction, i.e. when the coupling function / = 1, the system becomes uncoupled: the evolution 
of the fluid does not affect the distribution of the scalar field. This situation differs from the case of the presence of 
strong gravitational fields considered in Ref. 13] when even at / = 1 equations describing static configurations are 
not independent, but are coupled via gravity. 

In the case of the polytropic equation of state with polytropic exponent 7 

p = Kp 1 , (A.13) 

where K is some positive arbitrary constant, Eqs. (|A.9[) - (|A.12[> permit the introduction of similarity variables given 
by |H: 



v(r, t) = \fkt n - l u{x), p{r, t) = 



a{x) 



jL u3/2±3n-2 K \ ■ ) 

p(r,t) = —t 2 "- 4 [a(x)y, M(r,t) = - —m{x), <p(r, t) = -== <f>(x), 

AirG (3n - 2)G VAnG 

where k is some dimensional constant, and n - a dimensionless constant. Rewriting Eqs. (IA.9I) and (|A.10[) in terms 
of the above dimensionless similarity variables, we get an algebraic expression for m(x) in terms of a(x) and u(x) 

m(x) = f x 2 a(nx — u) . (A. 15) 

Next, Eq. (|A.11|) gives 

y-2^a . .du , nx-u 

7" -1 inx-u)— = -{n-l)u- }- -a. (A. 16) 

ax ax Sn — 2 

Differentiating now Eq. (|A.15|) with respect to x and taking into account (IA.10[) . we get 

. .da du x — u , .dhifdfj) , . _ s 

(nx — u)— a— = —2 a — a(nx — u) — — — . (A. 17) 

dx dx x do dx 
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Next, Eq. (jA~l~2|) takes the fo rm 



d 2 ^ 
dx 2 



2 deb 
x dx 



AirGk , 2 „ dV 



K 2 



-t 



^_/4(n-l) 7^. 
K 2 deb 



(A.18) 



It can be seen from Eqs. (|A.13[) and (|A.14[) that the above similarity variables permit the following relation: 

K = fc^TrG) 7 - 1 ^"^- 2 ). (A.19) 

Taking this expression into account, one can see from (|A.18j) that it is impossible to exclude the time variable t in 
front of dV/deb and df /deb simultaneously (as required by the self-similar character of the motion under consideration). 
That is why, bearing in mind that we look for solutions with the function / ^ 1, we will consider a massless scalar 
field with V{eb) = only. In this case the last term on the right-hand side of Eq. (|A.18|) is proportional to £ 4 ( 1 ~t). 
Hence the exponent is equal to zero when 7 = 1. Taking this into account and requiring that the parameter K from 
(|A.19[) be independent of time, i.e. that (n + 7 — 2) = 0, it is necessary to take n = 1. Then we finally have from 



f0 

dx 2 



2 deb 
x dx 



df 



Letting 7 = n = 1, Eqs. (|A.16[) and (|A. 17|i . in turn, can be rewritten in the following form: 

din / d- 



(x- 



a 

u)' 
1 



da 
dx 
du 
dx 

d 2 eb 2 deb 
dx 2 x dx 



fa- 2- 



(x — u)- 



(x 



2 din f deb 

fa{x-u) -J7--T 

x deb dx 



dx 
(x - 11), 



u), 



-a- 



df 



(A.20) 

(A.21) 
(A.22) 
(A.23) 



When the scalar field is absent, the system (IA.21[) - (|A.23[) describes the known problem of the self-similar motion of 



an isothermal fluid considered, e.g. in Refs. [15H18I 
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